Spectral functions and pseudogap in the t-J model 
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We calculate spectral functions within the t-J model as relevant to cuprates in the regime from low 
to optimum doping. On the basis of equations of motion for projected operators an effective spin- 
fermion coupling is derived. The self energy due to short-wavelength transverse spin fluctuations 
is shown to lead to a modifled selfconsistent Born approximation, which can explain strong asym- 
metry between hole and electron quasiparticles. The coupling to long-wavelength longitudinal spin 
fluctuations governs the low-frequency behavior and results in a pseudogap, which at low doping 
effectively truncates the Fermi surface, in particular near the (tt, 0) point in the Brillouin zone. 

PACS numbers: 7f .27.-)-a, 72.15.-v, 71.10.Fd 



In recent years underdoped cuprates are in the cen- 
ter of experimental and theoretical investigation, offer- 
ing possibly a clue to the understanding of anomalous 
normal state properties and the mechanism of supercon- 
ductivity in these compounds. Here we concentrate on 
some experimental facts revealing the nature of quasi- 
particles (QP) and the pseudogap. Several quantities, 
in particular the uniform susceptibility, the Hall con- 
stant, the specific heat, show the (large) pseudogap scale 
T* consistent with the angle resolved photoemission 
(ARPES) revealing a hump at ^ 100 meV observed in 
Bi2Sr2CaCu202+5 (BSCCO) near the (7r,0) momentum 
At the same time QP dispersing through the Fermi 
surface (FS) are resolved by ARPES in BSCCO only in 
parts of the large FS, in particular along the nodal (0, 0)- 
(tt, tt) direction, indicating that the rest of the large FS is 
either fully or effectively gaped. On approaching the op- 
timal doping the T* scale merges with Tc and the large FS 
becomes well defined. All these phenomena are naturally 
associated with the antiferromagnetic (AFM) short-range 
order (SRO) in cuprates, since the scale T* approaches 
the AFM exchange J in undopcd materials. This is less 
clear for the lower spin-gap scale Tgg (not the subject 
here), not found in all cuprates. 

While these facts allow for a qualitatively consistent 
scenario, the consensus on necessary prerequisites and 
moreover a satisfactory theoretical analysis are still miss- 
ing. The renormalization group studies of the Hubbard 
model 1^ indicate on the breakdown of the standard 
Fermi liquid and on the truncation of the FS. That such 
features also emerge from prototype models of correlated 
electrons has been confirmed in numerical studies of spec- 
tral functions in the Hubbard Q and in the t-J model 
which both show the appearance of the pseudogap 
at low doping. Some aspects of the pseudogap have been 
found in the spin-fermion models and studied phe- 
nomenologically in the Hubbard model § . 

Our aim is to capture these features within an analyti- 
cal treatment of a single band model. In the following we 
show that an effective spin-fermion model can be derived 



via equations of motion (EQM) and dividing the coupling 
into short and long-wavelength spin fluctuations an ap- 
proximation for the electron self energy can be found. 
We study the planar t-J model 



H — tijCj^Ci, 



+ J^(S,; 



• - ^n,nj), (1) 



where we take into account possible longer range hop- 
ping, i.e., besides = t for n.n. hopping also tij — t' 
for n.n.n. neighbors on a square lattice. We evaluate the 
single-particle propagator in this model explicitly tak- 
ing into account that fermionic operators are projected 
ones not allowing for the double occupancy of sites, e.g., 

4 = (1 -'^i,-s)cL- 

We use EQM directly for projected operators joj and 
represent them in variables appropriate for a paramag- 
netic metallic state with (S^) ~ and electron concen- 
tration (rii) = Ce = 1 — Ch, 

Pk,, H] = [(1 - Jc,]c^, + -L ^(2J^q + el_^) 



(2) 



where = — 4i7k — 4i'7^ is the bare band energy on 
a square lattice and 7k = (cosfca, -I- cos ky)/ 2, 7^ = 



cos kx cos ky . 



EQM for Cks can be used to construct approximations 
for the electron propagator G'(k, cj) [^,0, which can be 
represented as 



G(k,Lu) 



a 



UJ + ii~Ck- 5](k, w) ' 



(3) 



where the renormalization a = (l+c/i)/2 is a consequence 
of the projected basis, /i is the chemical potential, and Ck 
is the 'free' propagation term emerging from the EQM, 

Ck = -{{[cks,H],~clj+) -C = -47,1 i7k - ^mt'iL, (4) 
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where rjj = a + (So-Sj)/a and C is a constant. The 
central quantity for further consideration is the self en- 
ergy E(k, u) = ({Cks ; C+))Z'^/a, where tC^s = [5k., H] - 
CkCks, and only the 'irreducible' part of the correlation 
function should be taken into account in the evaluation 
of S. In finding an approximation for S we assume that 
we are dealing with the paramagnet with pronounced 
AFM SRO with the dominant wave vector Q = (tt, tt) 
and the AFM correlation length ^ > 1 with correspond- 
ing K = 1/^. We first note that EQM, Eq. (|), natu- 
rally indicate on an effective coupling between fermions 
and spin degrees. However, the role of short-range and 
longer-range spin fluctuations is quite different. 

In an undoped system AFM the spectral function of 
an added hole is quite well described within the self- 
consistent Born approximation (SCBA) where the 
strong hole-magnon coupling induced by the hopping t- 
term leads to a broad background representing the in- 
coherent hopping and a narrow QP dispersion governed 
predominantly by J. If we assume as a starting point an 
undoped Neel state as well as J < t EQM (||) directly 
reproduce the coupling equivalent to the holon-spin cou- 
pling within the SCBA. Note that within the Neel (Ising) 
state we have rji = 0, 772 = 1, — ±1/2, and a non- 
trivial coupling comes from transverse which can be 
represented via magnon excitations. Therefore by per- 
forming the decoupling of fermion and spin degrees in S 
and by using the identity Cjs — Cj_-sSj, we recover the 
standard SCBA equations. 

Our EQM formalism thus naturally leads to the SCBA 
in an undoped system. Still we note that in an isotropic 
AFM 7/1 ^ (but I771I < 1) which shghtly modificate 
SCBA results. Since the SCBA accounts well for proper- 
ties of a single QP in an AFM, we are not trying here to 
improve it. Within our approach we generalize the equa- 
tions for finite doping Ch > where we have electron-like 
QP above the Fermi energy {lj > 0). In 2D the AFM 
long-range order is absent due to T > and Ch > 0, still 
spin fluctuations are magnon-like, i.e., propagating and 
transverse to the local AFM SRO, with a dispersion ojq 
for q > K and q > k where q = q — Q. Hence the para- 
magnon contribution to the self energy can be written 
as 

Spm(k, W) = (Wq7k-q + V^-fl^f 

q,q>K, 

[G-(k^(l,Lj + ujq) + G+{k + q,uj-uj^)], (5) 



where (uq, Vq) = (1, -sign(7q))^(2 J ± Wq)/2tjq and 
refer to the Green's functions corresponding to electron 
{u! > 0) and hole {uj < 0) QP excitations, respectively. 
So far equations are written for T = 0, however in Epm 
the role of finite but low T > is not pronounced. Note 
that analogous to the SCBA t' does not enter directly 
the coupling but remains in the 'free' propagation term 
^k- Here we stress two features of our generalized SCBA: 



a) we are dealing with a strong coupling theory due to 
t > ujq hence a self consistent calculation of E is required, 

b) resulting spectral functions A(]s., to) are very asymmet- 
ric with respect to w = 0, since G"*" has less weight and 
consequently the scattering of electron QP is less pro- 
nounced. 

We are dealing with a paramagnet, therefore it is es- 
sential to consider also the coupling to longitudinal spin 
fluctuations. Note that the EQM (||) naturally introduce 
a coupling between fermion and spin operators which 
is isotropic in the spin space as appropriate in a para- 
magnetic state. In fact analogous form as in Eq. (|^) 
would emerge also from a spin-fermion Hamiltonian with 
the coupling parameter mtq = 2J7q -I- e^-q- ^"^^ 
fective Hamiltonian should be hermitian, i.e., the cou- 
pling should satisfy the condition Tfik.q = rhk-q.-q, 
therefore we use further on the symmetrized Wkq = 

2^7q+j(eLq + 4); 

Fermions and longitudinal spin fluctuations with q < k 
appear to be quite uncoupled, therefore we express the 
longitudinal contribution as in Refs. [p|jlC|], 



Eif(k,cj) 



1 \ ^ - 2 



aN 



dujiduj2 



9(^1,1^2) 



A0(k-q,u;i)x"(q,c^2) 



UJi — U>2 



(6) 



susceptibility, 
Enm)"^ and 



where x(q, cj) is the dynamical spin 
A°(k,u;) = -(Q:/7r)Im(w + - Ck - Sp„i)" 
5(tJi,W2) = th(wi/2T) -I- cth(tj2/2r). In Eif only the 
part corresponding to irreducible diagrams should en- 
ter, so there are restrictions on proper decoupling. We 
are mostly dealing with the situation with a pronounced 
AFM SRO where longitudinal spin fluctuations are slow, 
with a characteristic frequency Wk ^ J which is the 
case of a quasistatic x(q, w). Therefore we in Eq. (||) 
as the simplest approximation insert the unrenormal- 
ized A*'(k, a;), i.e., the spectral function without a self- 
consistent consideration of Eif but with Spm fully taken 
into account. Such an approximation has been intro- 
duced in the theory of a pseudogap in charge density 
wave systems p^], used also in related works analyzing 
the role of spin fluctuations ]T^, p4[ , and recently more 
extensively examined in Ref. |I5|. 

So far we do not have a corresponding theory for the 
spin response at c/i > and T > 0, so w) is assumed 
as a phenomenological input, bound by the sum rule 



cn). (7) 



The response should qualitatively correspond to a para- 
magnet close to the AFM instability, so we assume the 
form 



X"(q,t^) oc 



{q^ + n^)[bJ^ + uS^y 



(8) 
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where ^(w, T) oc lo would be appropriate for a nearly 
AFM Fermi liquid or an undoped AFM in 2D at 

any T > 0. On the other hand, in cuprates at inter- 
mediate doping more consistent with model results for 
T > seems to be the marginal Fermi liquid behavior 
with </>(w,T) a th(w/2T) |l6|,|l]§. 

Eqs. (|,|) for E — Spin + Sif represent the selfcon- 
sistent set of equations for G. Parameters k, 771, 772 are 
mainly dependent on Ch and are known from model cal- 
culations At T = and given Ch we determine fi 
such that the density of states J\f{uj) = (2/7V) A{k, uj) 
integrated for < reproduces Cg. At the same time FS 
is given by the relation Cicf + "^'O^f, 0) = /i. Full numer- 
ical analysis of selfconsistent equations will be presented 
elsewhere. Here we concentrate on some key aspects of 
the theory. One is that Epm allows for a meaningful be- 
havior in the limit Ch 0, which has been the deficiency 
of most phenomenological theories so far. In this limit 
our results for A(k, lj) at T ~ are essentially equivalent 
to results within the SCBA approach . For c^j > we 
have a finite contribution from electron QP A'^(k, a; > 0) 
and the corresponding QP density evolves as oc c/j. 

We choose further on parameters J = 0.3t, t' = —0.2t 
and K — 2y^. In Fig. 1 we present typical results for 
v4(k, oj) along the (0, 0) — (tt, tt) direction. Since at low 
doping r]2 0.9 main doping dependence arises from 
rii{ch), which varies from rji ~ —0.18 at Ch to rji ^ a 
for large doping. In calculation we use r]i close to values 
emerging from spin correlations found numerically |T7[ . 

As presented in Fig. 1 Epn^ leads to a strong damp- 
ing of the hole QP and quite incoherent momentum- 
independent spectrum A(k, w) for uj <^ — J which qual- 
itatively reproduces ARPES and numerical results 
Electron QP (at w > 0) are in general very different, i.e., 
with much weaker damping arising only from Spm. We 
should note that at given fi, Ce calculated from the den- 
sity of states does not in general coincide with the one 
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FIG. 1. ACkju) along the (0,0) — (7r,7r) direction for 
J = OM and t' = -Q.2t. (a) ch = 0.04, 771 = 0.02 and 
(b) Ch = 0.18, 771 = 0.2. 



evaluated from the FS volume, Cg = Vfs/Vq. Never- 
theless, apart from the fact that within the t-J model 
the validity of the Luttinger theorem is anyhow under 
question , in the regimes of large FS both quantities 
appear to be quite close. 

Results for a characteristic development of the FS with 
Ch are shown in Fig. 2. At Ch < c^^ ^ 0.08 solutions are 
consistent with a small pocket- like FS, whereby this be- 
havior is enhanced by <' < as realized in other model 
studies [|l^ . On increasing doping Ch > the FS rather 
abruptly changes from small into a large one. The small- 
ness of cjj has the origin in quite weak dispersion domi- 
nated by J and t' at c/i —> which is overshadowed by 
much larger (^k at moderate doping, where the FS is large 
and its shape is controlled by t' /t. Nevertheless, the rele- 
vance of obtained FS should be considered in connection 
with a coexistent pscudogap discussed further on. 

The position of the FS is mainly determined by Ck 
and Spin, while in this respect Sif is less crucial. Re- 
sults for Spjjj(k, 0) can be well parameterized in the 
form obtained within the SCBA for a single hole 
Similarly one can present also the full effective QP dis- 
persion, = Ck + S'(k, 0) — /i, and the QP residue 
Zk. The simplest approximation to discuss the pscu- 
dogap is the quasi-static approximation which is mean- 
ingful for Lu,^ <^ t. Assuming also k <C 1 simplifying 
,gx"(q, cj) ^ 7r(5(q — Q)5(u;)/4, as well as a single-pole 
form v4"(k, cj) = aZk(5(cj — e^) near the FS, we obtain 
from Eqs. (gi) 



G(k,w) 



where 



- ek-Q)(^ 

function 



^kQ 



(9) 



the gap function is given by A^q = 
■Z^k-Z^k-Q"^^q/4. From resulting branches it is ev- 
ident that a gap opens on the AFM zone boundary, so 
that the relevant pscudogap energy is A^*^ = | Ak^pM I ^ 
Zk|2J-4t'cos^/ca;|/2. For t' < the gap is largest at (tt, 0) 
point, as observed in experiments Since within the 
same k region the QP dispersion is also smallest the 
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FIG. 2. Fermi surface corresponding to results presented in 
Fig. 1. (a) Small FS, ch = 0.04 and (b) large FS, Ch = 0.18. 
Contour lines represent QP energy levels in increments Q.lt. 
The region with a developed pseudogap (w/A < 1) is 
line-shaded while the grey-shaded region represents the re- 
gion where the pseudogap is smeared out (w/A > 1). 
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effect is even more pronounced. If a; = is in the regime 
of the gap, then naturally we are dealing (in this approx- 
imation) with a truncated FS. For parameters as above 
we present in Fig. 2(b) the pseudogap region where the 
states and FS near the (tt, 0) momentum are strongly 
suppressed. 

Going beyond the quasi-static and k ^ treatment 
one can discuss also the QP states within the pseudogap. 
To study the general structure of the SF in this region it 
is enough to follow the development with e — e"^ crossing 
the pseudogap perpendicular to the AFM zone boundary. 
It is essential to take into account k > so that the 
averaging over q leads to an effective smearing of the 
delta-function ^^-q ^^1° ^ broader A(e,a;). So we have 
qualitatively to deal (at T = 0) with the self energy 



S"(e,tj)cx / x"iuj~ij')A{e,Lu')dij', 



(10) 



where the simplest assumption for A{e, uj) is a Lorentzian 
with the width w = v-^pK. Analogous equations have 
been already studied in Ref. |13| and lead to the pseudo- 
gap of the order of A = A^*^ , pronounced in QP spectra 
and clearly in J\f{uj). Results for the case with the gap 
centered at w = are shown in Fig. 3, where the deple- 
tion is most evident for w ^ A, while the pseudogap fills 
up for w > A. 
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FIG. 3. Density of states as a function of a; /A for 

= 0.2A. 

Looking at corresponding spectral functions directly, 
we notice that for a developed pseudogap with it; <C A (in 
Fig. 2(b) line-shaded) there are still QP crossing the FS, 
although with small ^ 1, while their velocity is not 
diminished. On the other hand if w > A the pseudogap is 
smeared out and consequently not effective, hence the FS 
is fully recovered. At intermediate doping this typically 
happens near the zone diagonal as shown in Fig. 2(b) 
(grey shaded). 

In conclusion, we have presented a theory for the spec- 
tral functions within the t-J model where the double- 
occupancy constraint is taken explicitly into account and 
used to derive an effective spin-fermion coupling. The 
coupling to transverse AFM paramagnons is strong, nev- 



ertheless it can be well treated within a generalized 
SCBA. On the other hand, the coupling to longitudinal 
AFM ffuctuations, mkq, is moderate near FS for low dop- 
ing and leads to a pseudogap, fully developed near the 
(tt, 0) point. The pseudogap is not in contradiction with 
the existence of a large FS, and should show up in in- 
tegrated photoemission and ARPES results as well as in 
the uniform susceptibility and in the specific heat. More 
elaborate analysis of proposed theory will be presented 
elsewhere. 
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